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We investigate the phase structure of two-color QCD at both real and imaginary chemical potentials {pf, 
performing lattice simulations and analyzing the data with the Polyakov-loop extended Nambu-Jona-Lasinio 
(PNJL) model. Lattice QCD simulations are done on an 8® x 4 lattice with the clover-improved two-flavor 
Wilson fermion action and the renormalization-group improved Iwasaki gauge action. We test the analytic 
continuation of physical quantities from imaginary p to real p by comparing lattice QCD results calculated at 
real p with the result of analytic function the coefficients of which are determined from lattice QCD results 
at imaginary p. We also test the validity of the PNJL model by comparing model results with lattice QCD 
ones. The PNJL model is good in the deconfinement region, but less accurate in the transition and confinement 
regions. This problem is improved by introducing the baryon degree of freedom to the model. It is also found 
that the vector-type four-quark interaction is necessary to explain lattice data on the quark number density. 
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I. INTRODUCTION 

Exploration of QCD phase diagram is one of the most im¬ 
portant subjects in not only nuclear and particle physics but 
also cosmology and astrophysics. However, due to the com¬ 
plexity of fermion determinant, the first principle calculation, 
i.e., lattice QCD (LQCD) simulations are quite difficult at 
high quark number chemical potential p. The QCD partition 
function Z at finite temperature T and finite p is expressed by 

Z = J DUdet[M{p)]e-^’’, (1) 

where (/r = 1, 2, 3,4) and Sq are the link variables and 
the pure gauge action, respectively, and M (p) is written as 

M{p) = + TO - PlA (2) 

with the covariant derivative and the quark mass to in the 
continuum limit. For later convenience, we regard p as a com¬ 
plex variable. It is easy to verify 

{det[M(p)]}* ** = det[M{-p*)]. (3) 


* 12634019@edu.cc.saga-u.ac.jp 
1 ysakai@riken.jp 

I tsaitou@kochi-u.ac.jp 
Siishii@email.phys.kyushu-u.ac.jp 

takahashi @ phy s .kyushu-u. ac .jp 

** kouji.kashiwa@yukawa.kyoto-u.ac.jp 

II kounoh@cc.saga-u.ac.jp 

11 nakamura @ riise .hiroshima-u .ac.jp 
Sl-yahiro@phys.kyushu-u.ac.jp 


Hence, the fermion determinant det[M{p)] is not real when 
p is real, and the importance sampling technique does not 
work in the Monte Carlo simulations there. This is the well- 
known sign problem. Several methods were proposed so far 
to resolve this problem; these are the reweighting method IH], 
the Taylor expansion method 0 ^ the analytic continuation 
from im aginary p to real p iRI- llffl , the complex Langevin 
method iflll - fT^ and the Lefschetz thimble theory llldl IlSll . 
However, these are still far from perfection. 

On the contrary, in two color QCD (QC 2 D), the lattice sim¬ 
ulations can be made at real and finite p, since the theory has 
no sign problem In fact, the following relation is 

obtained: 

det[M{p)] = det[{t2Cj5)~^M{p){t2C'y5)] 

= idet[Mip*)]r, (4) 

where t 2 and C = 7274 are the second Pauli matrix in color 
space and the charge conjugation matrix, respectively. Obvi¬ 
ously, det[M{p)] is real when p is real. Recently, Hands et 
al. analyzed the phase structure of QC 2 D in a wide range of 
real p by using two-flavor Wilson fermions QC 2 D 

can be also used to check the validity of methods proposed to 
resolve the sign problem. In fact, Cea et al. 11^ 1^ tested the 
validity of analytic continuation from imaginary p to real p 
with staggered fermions. 

Equation (|2l shows that det is real when p is pure 

imaginary, i.e., p = ipi = iOT for real variables pi and 6, so 
that LQCD simulations are feasible there. Observables at real 
p are extracted from those at imaginary p with analytic con¬ 
tinuation. In the analytic continuation, we must pay attention 
to the structure of phase diagram in the imaginary p region 
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where QCD has two characteristic properties, the Roberge- 
Weiss (RW) periodicity and the RW transition ll24l] . The QCD 
grand partition function has a periodicity of 27r/iVc in 9: 

Z«,) = z(9+^) ,5) 

for integer k and the number of color Nc. This periodicity was 
found by Roberge and Weiss and is then called the RW period¬ 
icity. Roberge and Weiss also showed that a first-order phase 
transition occurs at T > Trw and 9 = {2k -f l)7r/7Vc. This 
transition is named the RW transition, and T^w is slightly 
larger than the pseudo-critical temperature T^o of deconfine¬ 
ment transition at zero p. These features are remnants of Z 
symmetry in the pure gauge limit. These properties are con¬ 
firmed by LQCD simulations 

The RW periodicity does not mean that Zjsf^ symmetry is 
exact. Hence, there is no a priori reason that the order param¬ 
eter for ZjVc symmetry such as the Polyakov loop <P is zero 
in the confinement phase. In fact, in the case of Nc = 3, the 
Polyakov loop is always finite even in the confinement phase, 
when T is finite. However, the case of TVc = 2 is special 
In this case, the action and the boundary conditions are in¬ 
variant at /ii/T = {2k -f l)7r/2 under the CZ 2 transforma¬ 
tion composed of the Z 2 transformation and charge conjuga¬ 
tion C 1 I 25 I] . Due to this symmetry, the Polyakov loop becomes 
zero at low T when ^.i/T = {2k + 1)^/2. Pa ying attention 
to these characteristic features, Cea et al. 11^ 1^ analyzed 
the validity of analytic continuation in QC 2 D and found that 
lattice QC 2 D (LQC 2 D) data at real p can be described by a 
suitable analytic function, when the coefficients of analytic 
function are determined fromLQC 2 D data at imaginary /r. 

The results of LQCD at imaginary p are also useful 
to determine the parameters of effective models, such as 
the Polyakov-loop extended Nambu-Jona-Lasinio (PNJL) 
model 0261 - 1^ . Here we call this approach ’’imaginary chem¬ 
ical potential matching approach” 13311 . It is known that the 
PNJL model can reproduce the results of LQCD at imaginary 
/i, at least qualitatively, since the model has the RW periodic¬ 
ity and the RW transition ll^ [sH] . It was proposed | ^ | that 
the strength Gv of vector-type four-quark interaction 137. [^ . 
which is expected to be important for the physics of neutron 
stars, may be determined from LQCD data at imaginary /r; 
for the relation between neutron star properties and Gy, see 
Ref. and references therein. In Refs. ||4^ and ||4^ , in 
fact, Gy is determined with this prescription. The validity of 
such a determination of parameters in effective models can be 
checked in QC 2 D. 

In this paper, we study the phase structure of QC 2 D at 
both real and imaginary /r by performing simulations on 
an 8^ X 4 lattice with the renormalization-group improved 
Iwasaki gauge action lld^ l43l] and the clover-improved two- 
flavor Wilson fermion action ll4^ and analyzing the QC 2 D 
data with the PNJL model. We first test the analytic continu¬ 
ation from imaginary /r to real /r by comparing LQC 2 D data 
calculated at real /r with the result of analytic function the co¬ 
efficients of which are determined from L QC 2 D data at imagi¬ 
nary p. Such a test was tested in Refs. 11^|2^ with staggered 
fermions. Here the test is made with clover-improved Wilson 


fermions by assuming a polynomial series in the deconflne- 
ment phase and a Fourier series in the confinement phase. 

We second test the validity of the PNJL model by compar¬ 
ing LQC 2 D results with model ones. The PNJL model is good 
in the deconfinement region, but less accurate in the confine¬ 
ment region. This problem is improved by introducing the 
baryon degree of freedom to the model. It is also found that 
the vector-type four-quark interaction is necessary to explain 
QC 2 D data on the quark number density Uq. 

This paper is organized as follows. Section ^presents the 
lattice action and the parameter setting used in our LQC 2 D 
simulations. The definition of physical quantities is also pre¬ 
sented. In Sec. [nil the PNJL model is recapitulated. In 
Sec.HV] numerical results of LQC 2 D are shown and the ana¬ 
lytical continuation of physical quantity from imaginary /r to 
real p is tested. Comparison between LQC 2 D data and PNJL 
results are made in Sec. Section Iyj is devoted to a sum¬ 
mary. 


II. LATTICE SIMULATIONS 
A. Lattice action 

We use the renormalization-group-improved Iwasaki gauge 
action 5*0 ll4^l43| an d the clover-improved two-flavor Wilson 
quark action Sq 0441] defined by 

S = Sg + Sq, 

sg = -pY. (^0 Y ^ 

+ci E 

Sq — 'y ^ y ^ Qx^x,y<ii ) 

f=u,d x,y 

where q is the quark held, [3 = 4/p^, ci = —0.331, cq = 
1 — 8ci, and 

3 

Mx,y = Sxy — ~ li)Ux,i5x+l.v 

i=l 

+ (1 + li)Ul^i5x,y+i} 

-K{e^{l - 74)G^,4i5^+4,y 
+e '^(1 -I- 74)t^y,4<5a;^j^+4} 

(9) 

Here k is the hopping parameter, Fi^^, is the lattice held 
strength and - fl,,)/{8i) with the standard 

clover-shaped combination of gauge links. 

The coefficient Csw of clover term is determined by using 
a result obtained in a perturbative mean-field improved value 
Csw = P~^l^ ||431] with the plaquette P calculated in one-loop 
perturbation theory, P = 1 — 0.3154/3“^ for Nc = 2 l^ . 


( 6 ) 




(7) 

( 8 ) 
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B. Parameter setting for simulations 

We denote temporal and spatial lattice sizes as Nt and Ns, 
respectively. The Hybrid Monte-Carlo algorithm is used to 
generate full QC 2 D configurations with two-flavor dynami¬ 
cal quarks. The simulations are performed on a lattice of 
X Nt = 8^ X 4. The step size of molecular dynamics is 
6t = 0.02 and the step number of the dynamics is Nt = 50. 
The acceptance ratio is more than 95%. We generated 10,000 
trajectories and removed the first 5,000 trajectories as thermal- 
ization for all the parameter set. The relation of parameters k 
and /3 to the corresponding T/Tco is determined by finding 
the line of constant physics where the ratio of the pseudo¬ 
scalar (PS) meson mass nips to the vector meson mass at 
T = /r = 0 are invariant; see Table |T]for the relation. 


Ns 

Nt 

P 


r/Tco 

8 

4 

0.60000 

0.13782 

0.87783(585) 



0.64000 

0.13770 

0.94126(628) 



0.66000 

0.13751 

0.98577(657) 



0.68000 

0.13695 

1.02608(84) 



0.70000 

0.13677 

1.06282(709) 



0.72000 

0.13679 

1.11629(745) 



0.74000 

0.13567 

1.16526(77) 



0.76000 

0.13493 

1.20197(802) 



0.78000 

0.13443 

1.24298(829) 


TABLE I: Summary of simulation parameters. T^o is the pseudo- 
critical temperature at p = 0. In this parameter setting, the lattice 
spacing a is about 1.38 1.96 GeV”*^ and the ratio rUps/mv is 0.8 

at r = 0. 


C. Physical observables 


In this paper, we calculate the Polyakov loop (<?), the quark 
number density (ug) and the chiral condensate (cr). The quark 
number density is calculated by 


Td 

n„ = —— log z, 


V 


where V is the spatial volume, and a is by 


( 10 ) 


o- = {qq), 


( 11 ) 


where q is quark field and the (O) is the average value of 
physical quantity O. The chiral condensate suffers from the 
renormalization, and chiral symmetry is explicitly broken by 
Wilson fermions. This makes it difficult to deal with the ab¬ 
solute value of chiral condensate itself. We then consider a 
variation 


6a{T, fj,) = a{T, fj,) - a{T, 0). 
The Polyakov-loop operator is defined by 




Nt 


N, 


( 12 ) 


(13) 


with link variables Up, S SU(2). The average value ^ of L is 
related to the single static-quark free energy Fg as 

<p={L) (14) 

The Polyakov loop is an order parameter of the confine¬ 
ment/deconfinement transition if m is infinitely large. In fact, 
if Fg is finite (infinite), is finite (zero). The symmetry as¬ 
sociated with the confinement/deconfinement transition is Z 2 
symmetry under the transformation 

f/4(x,f) ^ Z2it)U4,{x,t), (15) 

where Z 2 is the element of Z 2 group that depends only on 
the temporal coordinate t. Pure gauge action is invariant un¬ 
der this transformation while L is not. Hence, F = (L) is 
an order parameter of Z 2 symmetry breaking. Effects of dy¬ 
namical quark break Z 2 symmetry explicitly and F is not a 
proper order parameter of confinement/deconfinement tran¬ 
sition. As mentioned in the previous section, however, at 
6 — (2k + 4 )tt/2 the system is symmetric under the CZ 2 
transformation 1^ . Hence, F becomes an order parameter of 
the combined symmetry there. 


III. PNJL MODEL 

Two-color QCD has Pauli-Giirsey symmetry in the limit of 
m = = 0 The PNJL Lagrangian of QC 2 D is so 

constructed as to have the symmetry and is given by 

C = q{iYDv — m)q 

+ G[(qqf + (qil^rqf + \q^ Ci-i5T2t2q\^] 

+ GsKqqf + (qilbfqf + \q^ Gi-i^T2t2q\^f 

- G,(qj''qf -U(F), (16) 

where q, m, U, n, G, Gs, Gv are the two-flavor quark fields, 
the current quark mass, the Pauli matrices in the color and the 
flavor spaces, the coupling constants of the scalar-type four- 
quark interaction, the scalar-type eight-quark interaction and 
the vector-type four-quark interaction, respectively. The po¬ 
tential W is a function of F. 

The mean-field approximation leads us to the thermody¬ 
namical potential 17 as ll4^ 

f? = -2Nf I ^^[E+ + E;+T(ln/-+ln/+)] 

-f U + U(F) (17) 

with 

= l + 2Fe-^^^ (18) 

U = G(a^ + A^) + 3Gs(<j‘^ + 

- G^n% (19) 

where Nf is the number of flavors, and A = | (q^ Ci^ 5 T 2 t 2 q) \ 

is the diquark condensate. The E'^ are defined by 

Ep = sgn(£'p ± jl)^{Ep±jiy F A^, 


( 20 ) 
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where Ep = with the effective quark mass M = 

m — 2Gcr — AGso'{a'^ + A"^), fi = p — 2Gvnq, A = —2GA, 
and sgn(£'p ± /I) is the sign function. When m = p = 0, f2 
becomes invariant under the rotation in a-A plane as a con¬ 
sequence of Pauli-Giirsey symmetry. Usually, G and G s ar e 
assumed to be constant. However, they may depend on (p 11^ . 
Here we consider ^-dependent G and Gg defined by 

G = Go(l-a<?2)^ Gg = G8,o(l-a^^), (21) 

where Gq, Gg^o a are constant parameters. In the 
Polyakov gauge, <l> is given by 

<?=i(e*^ + e-*^)=cos(,^), (22) 

for real number </). Following Ref. we take the Polyakov- 
loop effective potential of the form 


G(<?) 

T 


-b 


+ In (1 



(23) 


As will be mentioned in the next section, we determine these 
parameters to reproduce LQC 2 D data on rig at /3 = 0.75 and 
the pseudo-critical temperature Tco = 146 MeV at /r = 0. 

Because this model is nonrenormalizable and the first two 
terms in 17 are divergent. We then regularize them by intro¬ 
ducing a three-dimensional momentum cutoff as 



The mean fields X = a,A,nq,<P are determined from the 
stationary conditions 



(25) 


We found that A is always zero at the temperature and the 
density we treat in this paper. 

There are six parameters in the NIL sector of the PNJL 
model. The values of m, Gq, Gg 0 and A are chosen to repro¬ 
duce the vacuum property in LQC 2 D simulations. The values 
of Gv and a are chosen to reproduce LQC 2 D data on see 
Sec. |IV] Table HJ shows the parameter set used in this paper. 


mps [MeV] 

Mo [MeV] 

Go[GeV-^] 

G8,o[GeV-®] 

616 

354 

4.6 

60 

Gv/Go 

A[MeV] 

m[MeV] 

a 

0.15 

700 

110 

0.2 


TABLE II: Parameters of PNJL model. Mo is the effective quark 
mass at vacuum. 


IV. ANALYTICAL CONTINUATION OF PHYSICAL 
QUANTITIES 


A. Analytical continuation 

In this section, we show the numerical results of LQC 2 D 
simulations and perform the analytical continuation of physi¬ 


cal quantities from the region at imaginary p to the region at 
real /i, and finally examine the validity of analytical continua¬ 
tion. 

Figure [T] shows T dependence of ^ for several values of 
/t = fj,/T from iTr/2 to 1.2. For all the cases except = 
(i7r/2) = — (tt/2) , (p increases smoothly as T goes up. 

The deconfinement transition is thus crossover there. This 
property is the same as in QCD with three colors ifsHl . For 
each the pseudo-critical temperature (or /3c(A^)) 

is defined by the temperature where the susceptibility of 
becomes maximum. It is found from LQC 2 D simulations at 
/t = 0 that Tco = Tc(0) = 146 MeV. As for = (^7^/2)^ 
<P is almost zero below T = 1.12rc0 and increases rapidly 
above T = 1.14Tco. This indicates that the transition is the 
second order. As mentioned in Sec. HI at fi = i9 = in 12 the 
system is symmetric under the CZ 2 transformation, and as 
an order parameter of spontaneous breaking of the symmetry 
is zero below T^w and finite above Trw Therefore, T^w is 
located somewhere in a range of T = 1.12 ^ IAATcq. 


© 


P>.1.44 O 
fl ^.64 ^ 

/i^=0.16 + 

p^=0 X 
p ^=-0.394 ♦ 
p^=-1.58 ¥ 
p =(i;r/2)’n 




' ,5 i 


^ • 


00 EsSSlHllffllilSS B 


T/Tco 


Fig. 1: LQC 2 D results on T dependence of P for several values of 


Figure |2] shows the pseudo-critical line /3c (A^) in 
plane, and it is found that /3c(A^) decreases as increases. 
The value of Tc{fi^) at determined from 

/3c((i7r/2)^) is 1.16Tco and slightly larger than Trw estimated 
above. This inconsistency is understood as follows. In our 
calculations, /3 dependence of the Polyakov-loop susceptibil¬ 
ity is calculated every 0.01 for each value of and the data 
are fitted with a Gaussian function and the value of is 

evaluated as the point where the Gaussian function becomes 
maximum. The Gaussian fitting is, however, not valid in the 
vicinity of singular point (the endpoint of RW phase transi¬ 
tion). We then adopt the value 1.12 ^ l.ldTcO as Trw 
N ow we test the analytic continuation from imaginary /i 
to real p, for the cases of j3 = 0.75,0.70,0.65 and 0.60 that 
correspond to T/Tco = 1.18,1.06,0.96,0.88, respectively. 
The system is in the deconfinement (D) phase at /3 = 0.75, 
while it is in the confinement (C) phase at /3 = 0.60. At 
/3 = 0.70 (0.65), the system is in the C-phase when < 
— 1.15 (1.35) and in the D-phase when jj? > —1.15 (1.35). 
For each temperature, we then use different analytic functions 
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as explained below 
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Fig. 2: Pseudo-critical line of deconfinemnet transition in 
plane. In each of thin horizontal solid lines, /? is constant, and in each 
of thin vertical dotted lines is constant. At /? = 0.70, the left hand 
side of the left thin vertical dotted line belongs to to the C-Phase, and 
the right side does to the D-phase. At /3 = 0.65, the left hand side of 
the right thin vertical dotted line corresponds to the C-Phase, and the 
right side does to the D-phase. Note that Pa{Q) = 0.67. 


1. Trw < r(/3 = 0.75) 

At this temperature, due to the existence of RW transition, 
physical quantities cannot be described by any continuous pe¬ 
riodic function. Hence we use a polynomial series of the form 

A + (26) 

or 

A + Bfi'^ + Cfi\ (27) 

for /i-even quantities and a, where A, B, C are expansion 
coefficient. For a p-odd quantity Uq, we use 


Afi + 

(28) 

Afi + + Cfi^. 

(29) 


2. Tco < T{P = 0.70) < Trw 


At this temperature, the system is in the D-phase when 
> —1.15. We then use the same polynomial series 
as in the case of /3 = 0.75, but consider only the region 
— 1.15 < < Oasa fitting range. For (p and a, we use 

only the quadratic function ( |26] |. since the number of data we 
can use is small. 


3. T{l3 = 0.60, 0.65) < Tco 

At this temperature, the system is in the C-phase at imagi¬ 
nary and zero fi. Hence, it is expected that physical quantities 
can be well described by continuous periodic functions. Since 
‘P{0) is 0-even and has a periodicity of 27r in 0 = Im(/t), we 
use the following Fourier series 


Acos(0), 

(30) 

Acos(0) -F B cos(30), 

(31) 


for (p. Note that the terms of cos(20) and cos(40) as well 
as the constant term are excluded, since <P{t:/2 + 6') = 
— <P(7r/2 — O') for any 0'. The chiral condensate a{9) is a 
0-even and periodic function with a period tt. We then use the 
following Fourier series 

A + Bcos(20), (32) 

or 

A + Bcos(20) + Ccos(40), (33) 

for a. The quark number density nq{9) is a 0-odd and peri¬ 
odic function with a period tt. We therefore use the following 
Fourier series 

Asin(20), (34) 

or 

Asin(20)+ Bsin(40), (35) 

for Uq. Note that, in the case of /3 = 0.65, the system is in the 
D-phase when > 1.35. Hence, the Fourier series in which 
the coefficients are determined from LQCD data at imaginary 
fi and zero fi = 0 may not work there. 


4. Pseudo-critical line 

The pseudo-critical line is /i-even. We then use the 

polynomial series (l26l l and (l27l i. 


B. Quark number density 

First we consider the analytic continuation of Uq. Figured 
shows dependence of (uq/T)'^ for several values of T. 
The analytic continuation has errors coming from LQCD data 
at zero and imaginary fi. We then plot the upper and lower 
bounds of analytic continuation by a pair of same lines. The 
(uq/T)'^ are smooth at /i = 0, as expected. This is true for 
da and (p, as shown later. This guarantees that the analytic 
continuation from imaginary jl to real fi is possible. 

At /3 = 0.75 (T/Tco = 1.18), the system is in the D-phase 
and hence the polynomial series is used. The coefficients de¬ 
termined from LQCD data at imaginary /i are tabulated in Ta- 
ble|ni](a) of AppendixiAl The polynomial series up to /t^ well 
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reproduces LQC 2 D data in a wide range of 0 < < (1.2)2. 

Note that the analytic function deviates from LQC 2 D data 
near the first-order RW phase transition present at T > rnw 
and = -(7r/2)^- 

At A = 0.70 (T/Tco = 1.06), the system is in the C- 
phase when jj? < —1.15, while it is in the D-phase otherwise. 
Hence, we use only seven data in a range ofA^ = —1.15"^0 
to determine the coefficients of polynomial series. The coeffi¬ 
cients of the function are tabulated in TableHIIllbl of Appendix 
IaI LQC 2 D data calculated at real A lie between the upper and 
the lower bounds of the polynomial series up to A^ in ^ wide 
range of 0 < A^ < (l-2)^- 

At /3 = 0.65 {T/Tco = 0.96), the system is in the C-phase 
when fj? < 1.35, while it is in the D-phase otherwise. Hence, 
we use the Fourier series. The coefficients of the function are 
tabulated in Table HIJ (cl of Appendix lAl The analytic func¬ 
tions fail to reproduce LQC 2 D data calculated at real A when 
A^ > 0.4. The large deviation at large A^ may be originated in 
the fact that the system is in the D-phase there and the Fourier 
series may not be valid anymore. 

At A = 0.60 (T/Tco = 0.88), the system is in the C-phase. 
Therefore, we use the Fourier series. The coefficients of the 
function are tabulated in Table |III] (d) of Appendix The 
Fourier series up to the term sin (26) (sin (40)) are consistent 
with LQC 2 D data calculated at real A in a wide range 0 < 

A^ < 0.8(0 < A^ < (1-2)2). 

Comparing four cases of T/Tco with each other, one can see 
that the analytic continuation is reasonable at higher T/Tco 
where the system is always in the D-phase when A^ varies 
from — (7r/2)2 to (1.2)2 TjTco where the sys¬ 

tem is always in the C-phase when varies from —(7r/2)2 to 
(1.2)2. jsjg^j. T/Tco = 1, however, the system changes from 
the C-phase to the D-phase as A^ varies from —(7r/2)2 to a 
positive value. A simple analytic function cannot follow the 
complicated change properly. Therefore, the analytic contin¬ 
uation is reasonable except for the vicinity of deconfinement 
crossover. 


when [j? < 1.35, while it is in the D-phase otherwise. Hence, 
we use the Fourier series. The analytic functions is not consis¬ 
tent with LQC 2 D data calculated at real A when A^ > 0.4. As 
mentioned in the case of n,, this failure at large A may show 
that the system is in the D-phase there and the Fourier series 
becomes less reliable. At /3 = 0.60 (T/Tco = 0.88), the sys¬ 
tem is in the C-phase. Hence, we use the Fourier series. The 
analytic functions are consistent with LQC 2 D data calculated 
at real A in a range of 0 < A^ < 0.8 

D. Polyakov loop 

Figure |5] shows A^ dependence of <P at several values of T. 
Again, the upper and the lower bounds of analytic continua¬ 
tion are shown by a pair of same lines; see Table [HI] for the 
coefficients of analytic function. As for one can make the 
same discussion qualitatively as in Sec. II V Cl for Sa. 

E. Pseudo-critical line 

Figure |6] shows the transition line of deconfinement 
crossover in plane. The pseudo-critical /3c(A^) pP = 
0 is about 0.67. We consider the polynomial series the coef¬ 
ficients of which are obtained from LQC 2 D data at imaginary 
/j, and tabulated in Table |IV] of Appendix 0 The polynomial 
series up to well reproduces LQC 2 D data at 0 < A^ < 0.8, 
but deviates at A^ > 0.8. The polynomial series up to A^ is 
consistent with LQC 2 D data even at A^ > 0.8, but the differ¬ 
ence between the upper and lower bounds of analytic contin¬ 
uation is large. Therefore, we should consider that the ana¬ 
lytic continuation of the pseudo-critical line is reasonable at 

0 < A^ < 0.8. 


C. Chiral condensate 

Figure 0] shows p^ dependence of 5a for several values of 
T. Again, the upper and lower bounds of analytic continu¬ 
ation are shown by a pair of same lines; see Table HID for the 
coefficients of analytic function determined from LQC 2 D data 
at imaginary /r. As for Sa, one can made the same discussion 
as in Sec. II V BI for Uq, as shown below, although the analytic 
function taken in the C-phase is a cosine function. 

At /3 = 0.75 {T/Tco = 1.18), the polynomial series up to 
A^ well reproduces LQC 2 D data calculated at real A in a wide 
range of 0 < A^ < (1-2)^. At /3 = 0.70 (T/T^o = 1.06), 
the system is in the C-phase when p^ < —1.15, while it is 
in the D-phase otherwise. Hence, we can use only three data 
in a range of A^ = —1.15 ^ 0 to determine the coefficients 
of analytic function and then use the quadratic function only. 
The function is consistent with LQC 2 D calculated at real p in 
a wide range of 0 < A^ < (l-2)^- 

At /3 = 0.65 {T/Tco = 0.96), the system is in the C-phase 
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/?=0.75 



8=0.70 



Fig. 3: dependence of Uq at (a) /? = 0.75 {T/Tco = 1.18), (b) 

P = 0.70 {T/Tco = 1.06), (c) P = 0.65 (T/Tco = 0.96) and (d) 
P = 0.60 (T/Tco = 0.88). The dots with error bars are the results 
of LQC 2 D data. The solid and dashed lines represent the results of 
analytic continuation in which two types of analytic functions are 
taken as shown by legends. The upper and lower bounds of analytic 
continuation is shown by a pair of same lines. Characters C and D 
denote confinement and deconfinement phases, respectively. 


^=0.75 



^=0.70 



ji =0.65 



^=0.60 



Fig. 4: -dependence of a at (a) /3 = 0.75 {T/Tco = 1.18) (b) 

P = 0.70 {T/Tco = 1.06) (c) P = 0.65 (T/T^o = 0.96) (d) 
P = 0.60 (T/Tco = 0.88). For the definition of dots, lines and 
characters, see the caption of Fig.[3 
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P =0.75 



=0.70 



^=0.60 



Pseudo-critical line 



Fig. 6: Pseudo-critical line of deconfinement crossover in 
plane. Note that /3c(0) = 0.67. For the definition of dots, lines and 
characters, see the caption of Fig. [S] 


Fig. 5: -dependence of $ at (a) /3 = 0.75 {T/Tco = 1.18), (b) 

13 = 0.70 (T/Tco = 1.06), (c) /3 = 0.65 {T/Tco = 0.96) and 
(d) P = 0.60 (T/Tco = 0.88). For the definition of dots, lines and 
characters, see the caption of Fig.[^ 
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V. COMPARISON OF PNJL MODEL WITH LQC 2 D DATA 

A. Parameter setting 

In this section, we compare results of the PNJL model with 
LQC 2 D data to test the validity of the model. For this purpose, 
we first fix the parameters of the model. For the NIL sector, 
the parameters have already been determined in Sec. |III] We 
then fix the remaining parameters a, b, a and Gv here. 

Figure|2]shows T dependence of Ug divided by its Stephan- 
Boltzmann (SB) limit nsB for several values of from 
— (7r/2)^ to (1.2)^. LQC 2 D results include a lattice artifact 
due to finite volume and spacing. The artifact is expected 
to be reduced in Ug/nsB- For all the values of (1^, the ratio 
Ug/nsB increases as T increases. 



Cr 

C 


fl ^=0.16 + 

fl ^=-0.394 X 
f'=-1.S8 X 

X 

♦ 


m ® & * 


A T -t- 


@ 

^ A 


T/Tco 


Fig. 7: T dependence of rig/nsB for several values of p,^. Six cases 
of , —1.58, —0.394, 0.16, 0.64,1.44 are shown from 

the bottom. 

Now we determine the parameters a, b, a and Gv from Ug 
at the highest T in the present analyses, i.e., at /3 = 0.75 
(T/Tco = 1.18). One reason is that the PNJL model is essen¬ 
tially a model for quark dynamics and it may work better at 
higher T than at lower T. Another reason is that Ug does not 
need the renormalization and is sensitive to the value of Gv. 

Figure m shows dependence of Ug/nsB at /3 = 0.75 at 
imaginary /r. As decreases, the chiral symmetry breaking 
becomes stronger and the effective quark mass M becomes 
larger 1^ . so that Ug/nsB decreases. We searched the pa¬ 
rameters a, 6, a and Gv so as to reproduce both the result 
of Fig. [Hand Tco = 146 MeV. The parameters thus obtained 
are shown in Table [III of Sec.|III| It is interesting that the vec¬ 
tor interaction is needed to reproduce LQC 2 D results data at 
imaginary /r. In fact, the model with Gv = 0.15Go (solid line) 
yields better agreement with LQC 2 D data than the model with 
Gv = 0 (dotted line). This method may work as a way of de- 
termining the vector coupling also in realistic QCD with three 
colors SSlill]. 

T dependence of <P is shown in Fig. |9]for three cases of 
p'^ = —(7r/2)^, 0,1.44. Since the renormalization is needed 


/3=0.75 



Fig. 8: p^ dependence of Ug/nsB at /3 = 0.75 at imaginary The 
dots with error bars are the results of LQC 2 D data. The solid (dotted) 
line is the result of PNJL model with Gv/Go = 0.15 (0). 


for <P, we multiply the PNJL results by a factor 0.304 to repro¬ 
duce LQC 2 D results at /t = 0 and T = The renormal¬ 
ized PNJL results (solid lines) reproduce LQC 2 D data qualita¬ 
tively, except for the vicinity of the first-order RW phase tran¬ 
sition at T > Tftw ~ l-13Tco and p^ = — (7r/2)^. The value 
of Trw is 178 MeV in the PNJL model, but 163 - 166MeV 
in LQC 2 D data. In the present model, it is quite difficult to 
reproduce LQC 2 D values of Tco and Trw simultaneously. A 
fine tuning of the Polyakov potential U{<P) may be necessary. 


0.5 I-^^-T-T- 

;i^=1.44 + PNJL - Ji^=l-44 

=0 V PNJL - P' =0 

/i'= (i?r/2)^ O p[,jj|_ . /i^=(in/2)^ 



-0.1 I-^^^-'-1 

0.8 0.9 1 1.1 1.2 1.3 

T/Tco 


Fig. 9: T dependence of # for three values of p^. The dots represent 
LQC 2 D data. The solid (dotted, dashed) line represents the PNJL 
results at = 0 (— (7r/2)^), 1.44). The PNJL results are multiplied 
hy the normalization factor 0.304. 


B. Quark number density 

Figure [TOj shows dependence of rig/nsB at ,8 = 
0.75 {T/Tea = 1-18). As mentioned in the previous subsec- 
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tion, we have used LQC 2 D on rig/nsB at /? = 0.75 and imag¬ 
inary /i to determine the parameter set of the PNJL model. 
The parameter set thus determined well reproduces LQC 2 D 
data even at real p. This ensures the assumption that the PNJL 
model is valid at high T. Also note that the PNJL model with 
Gv = 0 fails to reproduce LQC 2 D data at real /r, while the 
PNJL model with Gv agrees with the LQC 2 D data even at 
real p. The imaginary chemical potential matching approach 
is thus a promising method. 


P=0J5 



Fig. 10: -dependence of Uq/nss for /3 = 0.75 {T/Tco = 1.18). 

The dots with error bars are the results of LQC 2 D data, while the 

solid (dotted) line corresponds to the result of PNJL model with 

Gv/Go = 0.15 (0). Note that Uq/nsB at p = 0 is defined by 

lim Uq/nsB- 
p,^o 


Figure [TT] shows the same as Fig. m but for P = 
0.70 {T/Tco = 1.06). On the left (right) side of the vertical 
thin dotted line, the system is in the C-phase (D-phase). The 
PNJL result (solid line) underestimates LQC 2 D data sizably 
in the C-phase. To improve this, we consider free “baryons” 
and assume that the baryons have the same degrees of free¬ 
dom as PS mesons, and add a contribution of free baryon gas 
to the PNJL results. Note that, in QC 2 D, the baryon is a boson 
consisting of two quarks. According to Ref. jl^. the scalar 
baryon with the same degree of freedom as the PS meson has 
almost the same mass as the PS meson. Flence, we use the 
baryon mass me = rups = 616MeV. In this way, the baryon 
contribution to Uq is given by 


ng,B = 2g 


(fp 


1 

g/3(£^B-2Ai) 


1 


1 

g/3(£;B-l-2/i) _ I ’ 

(36) 


where Eb = \/p = 3 is the degree of freedom and 
the factor 2 in front of g comes from the fact that the baryon is 
composed of two quarks. This modification improves agree¬ 
ment with LQC 2 D particularly in the C-phase, but not in the 
D-phase. This implies that baryons disappear in the D-phase 
at least partially. 

Figure [12] shows the same as Fig. M but for (3 — 
0.65 {T/Tco = 0.96). The system is in the C-phase on the 



Fig. 11: -dependence of Uq/nsB for /J = 0.70 (T/Tco = 1.06). 
The dots with error bars are the results of LQC 2 D data, while the 
solid (dashed) line corresponds to the result of the PNJL model (the 
PNJL-l-baryon model). On the left (right) side of the thin dotted line, 
the system is in the C-phase (D-phase). 


left side of the thin dotted line, but in the D-phase on the 
right side. Again, the PNJL model undershoots LQC 2 D in the 
C-phase, but the PNJL-tbaryon model almost reproduces the 
LQC 2 D data in the C-phase, although the latter model over¬ 
shoots LQC 2 D data in the D-phase. Thus, the baryon may 
disappear in the D-phase. 

More precisely, the PNJL-tbaryon model overestimates 
LQC 2 D data also in the C-phase near the thin dotted line. This 
fact may imply that a repulsive force works between baryons 
there. It is well known that such a repulsive force suppresses 
the baryon number density in realistic nuclear matter. 



Fig. 12: -dependence of Uq/nsB for (5 = 0.65 (T/Tco = 0.96). 
For the definition of lines, see the caption of Fig. m 

Figure [T3] shows the same as Fig. [TOl but for /3 = 
0.60 {T/Tco = 0.88). The system is in the C-phase and the 
PNJL model largely underestimates LQC 2 D data, but it is im¬ 
proved by the PNJL-tbaryon model. We can therefore con¬ 
clude that baryon effects play an important contribution to Uq, 
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whenever the system is in the C-phase. 



Figures [15] and [T6| show -dependence of 5cr at ^ = 
0.70 (T/T^o = 1-06) and /3 = 0.65 {T/T^o = 0.96), re¬ 
spectively. On the left side of the thin vertical dotted line the 
system is in the C-phase, while it is in the D-phase on the right 
side. The PNJL model (solid line) is consistent with LQC 2 D 
data in the D-phase, but not in the C-phase. In order to im¬ 
prove this disagreement in the C-phase, we add baryon effects 
to the PNJL model again: 




dniB f dPp me 

dm ^ J (27r)3 Eb 
1 

g/3(t5B+2M) — 1 ’ 


1 

^l3(Es-'2.E) 


1 


(37) 


where we assume 


Fig. 13: -dependence of rig/nsB for /3 = 0.60 (T/Tco = 0.88). 
For the definition of lines, see the caption of Fig. M 


C. Chiral condensate 

As mentioned in Sec. Ill Cl the renormalization is necessary 
for the chiral condensate. The PNJL results are then simply 
multiplied by a normalization factor 1.92 so that the results 
can reproduce LQC 2 D data at /3 = 0.75 (T/Tco = 1.18) and 
= {in12)'^. This choice of normalization is natural, since 
the PNJL model is a quark model and is expected to be more 
reliable at higher T. 

Figure M shows dependence of (5cr at /3 = 

0.75 {T/Tco = 1.18). At this temperature, the system is in 
the D-phase, and the PNJL model well reproduces LQC 2 D 
data even at real /r. 



Fig. 14: -dependence of SajT^ at /3 = 0.75 (T/Tco = 1.18). 
The dots with error bars are the results of LQC 2 D data, while the 
solid line corresponds to the result of PNJL model. The PNJL result 
is multiplied by a normalization factor 1.92. 


dmB _ 2 
dm 


(38) 


since a naive constituent quark model gives this value. As 
shown in Figs [15] and [16] the PNJL-Hbaryon model (dashed 
line) is more consistent with LQC 2 D data than the PNJL 
model (solid line) in the C-phase, but less consistent in the 
D-phase. This means that baryon effects are significant only 
in the C-phase. 



Fig. 15: -dependence of da at /? = 0.70 {T/Tco = 1.06). The 
dots with error bars are the results of LQC 2 D data, while the solid 
(dashed) line corresponds to the results of the PNJL (PNJL+baryon) 
model. The model results are multiplied by a normalization factor 
1.92. The system is in the C-phase on the left side of the thin vertical 
dotted line, but in the D-phase on the right side. 


Figure \n\ shows -dependence of Sa at /3 = 0.60 
(T/Tco = 0.88). At this temperature, the system is in the 
C-phase. The PNJL-Hbaryon model (dashed line) yields better 
agreement with LQC 2 D than the PNJL model (solid line) in 
the whole region of note that 6a is always zero at /t = 0 
by the definition ([T2l i. Thus, baryon effects are important in 
the C-phase not only for Ug/nsB but also for Sa. 
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/3=0.65 




/3=0.75 



Fig. 16: -dependence of <5cr at /I = 0.65 (T /Tco = 0.96). For the 
definition of lines, see the caption of Fig. [B] 


/3=0.60 



Fig. 17: -dependence oi Sa/T^ at /3 = 0.60 (T/Tco = 0.88). 
For the definition of lines, see the caption of Fig. [B] 


Fig. 18: -dependence of <15 at /I = 0.75 (T/Tco = 1.18). The dots 
with error bars are the results of the LQC 2 D data, while the solid 
line corresponds to the result of PNJL model. The PNJL result is 
multiplied by the normalization factor 0.304. 



D. Polyakov loop 

FigurefTSlshows dependence of at ^ = 0.75. As men¬ 
tioned in Sec. IV Al the PNJL result is multiplied by a normal¬ 
ization factor 0.304 for <l>. At this temperature, the system is 
in the D-phase. The PNJL model (solid line) well reproduces 
LQC 2 D data except for the vicinity of the first-order RW tran¬ 
sition. 

Figures [T^ and l20l show dependence of at /3 = 0.70 
and /3 = 0.65, respectively. The PNJL model deviates from 
LQC 2 D data to some extent in the C-phase except for the point 
at = (i-|) where Z 2 symmetry restricts to zero. 

Figures I 2 T] shows -dependence of ^ at = 0.60. At this 
temperamre, the system is in the C-phase. The PNJL model 
deviates from LQC 2 D data to some extent in the whole region 
except for the point at = (i-KjZf' where symmetry 
restricts to zero. 

Throughout all the analyses for we can say that the PNJL 
model cannot reproduce LQC 2 D data properly in the vicinity 
of RW transition and at lower T. A fine tuning of 14 pP) may 


Fig. 19: -dependence ofPatP = 0.70 (T/Tco = 1.06). The dots 
with error bars are the results of the LQC 2 D data, while the solid 
line is the result of PNJL model. The PNJL result is multiplied by 
the normalization factor 0.304. The system is in the C-phase on the 
left side of the thin vertical dotted line, while it is in the D-phase on 
the right side. 


be necessary to solve this problem. 


E. Phase diagram 

Figure |22] shows the phase diagram in p^-T plane. The 
PNJL model (solid line) reproduces LQC 2 D data well in a 
range of p^ = — 2 ^ 1, but it overshoots LQC 2 D data to some 
extent near the RW transition line and undershoots LQC 2 D 
data in the large p"^ region of ~ (1-2)^. It is an interesting 
question whether the PNJL model can reproduce LQC 2 D data 
in the large p"^ region as soon as the model is improved to 
reproduce LQC 2 D data near the RW transition line. 
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(5 =0.65 



Phase diagram 



Fig. 20: //^-dependence of at /3 = 0.65 {T/Tco = 0.96). For the 
definition of lines, see the caption of Fig. EH 



Fig. 21: //^-dependence of ^ at /3 = 0.60 {T/Tco = 0.88). For the 
definition of lines, see the caption of Fig. 1181 


VI. SUMMARY 

We studied the phase structure of QC 2 D at both real 
and imaginary // by using an 8^ x 4 lattice with the 
renormalization-group improved Iwasaki gauge action iiS 
I 43 I] and the clover-improved two-flavor Wilson fermion ac¬ 
tion lld^ . The Polyakov loop, the chiral condensate and the 
quark number density were calculated at 0.86 < T/Tco < 
1.18 and —(7r/2)^ < ^ These quantities are 

smooth at A = 0, as expected. This guarantees that the an¬ 
alytic continuation of physical quantities from imaginary A to 
real A is possible. 

Accuracy of the analytic continuation was tested in Refs. 
ll^l^ with staggered fermions. In this paper we have made 
similar analyses with clover-improved Wilson fermions by as¬ 
suming a polynomial series of A in the deconfinement phase 
and a Fourier series in the confinement phase, where coeffi¬ 
cients of the series were determined at imaginary //. As for the 
quark number density at T jT^a = 1.18 corresponding to the 


Fig. 22: Phase diagram in jl^-T plane. The dots with error bars 
represent the Pseudo-critical line of LQC 2 D, while the solid line cor¬ 
responds to the phase diagram of PNJL model. The vertical thick 
solid line is the RW transition line determined with the PNJL model. 

deconfinement phase, the polynomial series up to A^ well re¬ 
produces LQC 2 D results in a wide range of 0 < A^ < ( 1 - 2 )^. 
At T/Tea = 0-88 corresponding to the confinement phase, 
the results of the lowest-order Fourier series sin(20) are con¬ 
sistent with LQC 2 D results in a range of 0 < A^ < 0 . 8 . At 
T/Tco = 0.96 near the deconfinement transition, it is good 
only in 0 < A^ < 0-4- The analytic continuation is thus 
useful in the deconfinement and confinement regions, but less 
accurate in the transition region near T/T^o — 1 where the 
deconfinement crossover takes place somewhere in a range of 
— ( 7 r/ 2 )^ < fP < (1.2)^ as jP increases with T fixed. This 
is true for other quantities such as the Polyakov loop and the 
chiral condensate. 

We have tested the validity of the PNJL model by compar¬ 
ing model results with LQC 2 D ones, where the model param¬ 
eters are fitted to the quark number density at T/Tco = 1.18 
and imaginary //. As for the transition line of deconfienment 
crossover, the model result agrees with LQC 2 D one. More 
precisely, the agreement is not perfect in the vicinity of RW 
transition line and the large-A^ region of (A)^ ~ (l-2)^- It 
is interesting whether the PNJL model can reproduce LQC 2 D 
data in the large (P region as soon as the model is improved 
to reproduce LQC 2 D data near the RW transition line. A 
possible candidate of the improvement is a fine tuning of the 
Poly akov-loop potential U{T>). 

In the deconfinement region of T/Tcq = 1.18, the PNJL 
model yields good agreement with LQC 2 D data at both real 
and imaginary // for the quark number density, the chiral con¬ 
densate and the Polyakov loop. The agreement particularly 
at real // indicates that the PNJL model is reliable in the de¬ 
confinement region. In the transition region of T/Tco ~ 1, 
the agreement of the PNJL model with LQC 2 D data is not 
perfect. As for the quark number density and the chiral con¬ 
densate, however, the deviation can be reduced in the con¬ 
finement area appearing at smaller jP by introducing baryon 
degree of freedom to the PNJL model. In the deconfinement 
area appearing at larger jP, on the contrary, the model overes- 
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timates LQC 2 D results if the baryon contribution is taken into 
account. This means that baryons disappear at least partially 
in the deconfinement area. Also in the confinement region of 
T jTca = 0.88, the baryon degree of freedom is important. As 
for the Polyakov loop, the disagreement between PNJL and 
LQC 2 D results in the confinement area cannot be solved by 
the baryon contribution. Of course, this comes from the fact 
that the Polyakov-loop potential is not changed by the baryon 
contribution in the present framework. The improvement of 
PNJL model along this line is interesting. 

The present analysis also shows that the vector-type four- 
quark interaction is necessary to explain LQC 2 D data on the 
quark number density. This fact indicates that also in the 
realistic case of three colors the strength of vector-type in¬ 
teraction can be determined from LQCD data at imaginary 
/i 1^Iddll . Although the lattice we used is quite coarse 
and the number of LQC 2 D data is limited, the present results 
surely show that the imaginary-p, matching approach is a 
promising approach to thermodynamics of QCD at finite p. 

In this paper, we consider only the p-region where the di¬ 
quark condensate is zero. Quantitative check of our effective 
model in the region lfl9l - l^ is also interesting. 
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Appendix A: Coefficients of fitting functions for analytic 
continuation 

We present the coefficients of analytic functions determined 
from LQC 2 D data at imaginary p and the x^/d.o.f for the 
fitting. 


(a) 

^ = 0.75 





observable 

function 

A 

B 

C 

xVd.o.f 

Im{nq) 

Eq.|2D 

3.61014(1466) 

0.84703(2265) 


2.914 

Im{nq) 

Eq.og 

3.57904(1989) 

0.70990(6347) 

-0.08459(3657) 

2.566 

Sa 


-0.07116(5561) 

-1.07870(5010) 


2.337 

5(7 

Eq.HT) 

0.01304(8702) 

-0.7336(0190) 

0.1571(8106) 

1.448 

# 

Eq.|26l 

0.26153(70) 

0.04379(12) 


13.25 


Eq.07) 

0.25994(73) 

0.02410(310) 

-0.01331(191) 

4.534 

(b) 

13 = 0.70 





observable 

function 

A 

B 

C 

xVd.o.f 

Im{nq) 

Eq.|2D 

3.32812(2089) 

1.07313(5618) 


1.077 

Im{nq) 

Eq.|22) 

3.34986(2818) 

1.2644(1755) 

0.24168(21017) 

1.016 

5(7 


0.00868(13029) 

-1.39546(25909) 


0.018 

# 

Eq.HT) 

0.21241(112) 

0.04621(279) 


5.590 

(c) 

0 = 0.65 





observable 

function 

A 

B 

C 

xVd.o.f 

Im{nq) 

Eq.O 

1.31566(759) 



23.80 

Im{nq) 

Eq-OD 

1.44376(1216) 

-0.11951(886) 


4.075 

5(7 

Eq.OH 

1.1107(5784) 

-1.33556(6933) 


2.155 

5(7 

Eq.03 

1.05958(6059) 

-1.27686(7235) 

0.19765(6969) 

0.1933 


Ea.im 

0.15644(85) 



7.920 

# 

Eg.im 

0.15375(94) 

-0.00461(68) 


0.442 

(d) 

= 0.60 





observable 

function 

A 

B 

C 

xVd.o.f 

Im{nq) 

Eq-OU 

0.93599(80) 



1.189 

Im{nq) 

Eq.OD 

0.94541(838) 

-0.01438(74) 


0.877 

5(7 

Eq.oa 

0.69931(3257) 

-0.71395(3973) 


0.2853 

5(7 

Eq.03 

0.69633(3414) 

-0.71059(4087) 

0.01373(3913) 

0.3394 

# 

Eq.oni 

0.10258(38) 



3.826 


Eq.BTI 

0.10404(52) 

-0.00195(46) 


0.4438 


TABLE III: Coefficients of analytic functions and x^/d.o.f for quark 
number density, chiral condensate and the Polyakov loop at (a) P = 
0.75, (b) P = 0.70, (c) P = 0.65, and (d) P = 0.60. The coefficients 
are determined from LQC 2 D data at imaginary p. The fitting range 
is p^ = —1.15 ~ 0 for (b) and p^ = —(7r/2)^ ~ 0 for the other 
cases. 


function ABC x^/do.f 
Eq.(|26j 0.66802(204) -0.02868(167) 0.109 

Eq.i|^ 0.66871(257) -0.02637(557) 0.00113(260) 0.089 


TABLE IV: Coefficients of analytic functions and x^/d.o.f for the 
pseudo-critical line. The coefficients are determined from LQC 2 D 
data at imaginary p. The fitting range is p^ = —( 77 / 2 )^ ~ 0 for all 
the cases. 
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